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We first reformulate and expand with several novel findings some of the basic results in the 
integrability of Abel equations. Next, these results are applied to Vein’s Abel equation whose 
solutions are expressed in terms of the third order hyperbolic functions and a phase space analysis 
of the corresponding nonlinear oscillator is also provided. 
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I. INTRODUCTION 


The first order Abel nonlinear equation of the second kind has the form 

, ,dw 2 s 

(iu + s)-- \-p + q 3 w + q 2 w + rw = 0 , (1) 

ax 

where p , q \, q\, r, and s are all some functions of x. The case r = 0 has been introduced more than two hundred 
years ago by Abel [lj . Many solvable equations of this type are collected in [1] and other ones can be found in more 
recent works M.' The transformation 1 /y = w + s converts this equation in the form 

= r+(q 2 - 3 rs)y + (qi - - 2q 2 s + 3 rs 2 )y 2 + (p - qis + q 2 s 2 - rs 3 )y 3 , (2) 

ax ax 

which is Abel’s equation of the first kind. We see that Abel’s original equation of the second kind is actually a 
homogeneous case of that of the first kind. 

One can also write (O in the canonical form 


-T- = fo( x ) + fi(x)y + f 2 (x)y 2 + f 3 (x)y 3 = F(x,y) , (3) 

ax 

where only the case s = 0 leads to the simple identifications: /o = r, /i = < 72 , fi = <7i, and f 3 = p. 

The integrability features of Abel’s equation in the canonical form are extremely important because of its connection 
with nonlinear second order differential equations that phenomenologically describe a wide class of nonlinear oscillators 
which in this way can be treated analytically. In particular, Vein [10] introduced an Abel equation with /o = 0 and 
rational forms of / 1 , f 2l and f 3 whose solutions are expressed in terms of third order hyperbolic functions. However, 
to the best of our knowledge, the paper of Vein went unnoticed for many years and only recently Yamaleev [Til, E3 
provided a generalization in the framework of third order multicomplex algebra. Needless to say, the corresponding 
Vein’s nonlinear oscillator has not been studied in the literature. This was the main motivation for writing this paper, 
which is organized as follows. In section [TT1 we present a lemma that encodes the connection of Abel equation with 
nonlinear oscillator equations and some related simple results. In section nm we reformulate several integrability 
results for Abel’s equation with the purpose to apply them to Vein’s Abel equation, which is the subject of section 
IIVI The dynamical systems analysis of Vein’s oscillator is developed in sectionfVj and we end up with the conclusions. 
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II. CONNECTION WITH SECOND ORDER NONLINEAR ODES 


The known importance of Abel’s equation in its canonical form ([3]) stems from the fact that its integrability leads 
to closed form solutions to a general nonlinear ODE of the form 

ZCC + h( x ) x c + f 3 (x) + fi(x)x 2 + fo(x)x 3 = 0 , (4) 

where the variables x(0 and y(x(0) are some parametric solutions that depend on a generalized coordinate This 
can be expressed by the following lemma Q 

Lemma 1. Solutions to a general second-order ODE of type fi} may be obtained via the solutions to 
Abel’s equation © and vice versa using the following relationship 

^ = «(*( 0 ) ■ ( 5 ) 


Proof. To show the equivalence, one just needs the chain rule 


d 2 x dv dx dv 

dC\ 2 dx d£ dx 


( 6 ) 


which turns © into the Abel equation of the second kind in canonical form 


V ^ + f 3 {x) + f 2 (x)v + fi(x)v 2 + fo{x)v 3 = 0. 


( 7 ) 


Via the inverse transformation 


z(z(C)) 


1 

2 / 0 ( 0 ) ' 


( 8 ) 


of the dependent variable, equation © becomes ©■ Moreover, the linear term in © can always be eliminated via 
the transformation 


y{x) = z{x)et h{x)d: 


( 9 ) 


which gives 


where 


dz 

dx 


h 0 (x) + h 2 {x)z 2 + h 3 (x)z 3 = H(x, z) , 


( 10 ) 


h 0 (x) = fo(x)e-f^ dx 

h 2 (x) = f 2 (x)el (11) 

h 3 (x) = f 3 (x)e 2 f f 'M dx . 


The case where /o = 0 —> r = 0 is seen from © to be the one actually considered by Abel, for in this case the 
reduced Abel equation 


, . dw n 

(w + s)— - \-p + qiiv + q 2 w =0 

dx 


( 12 ) 


can be always put in the form 


dz 

dx 


h 2 {x)z 2 + h 3 (x)z 3 


(13) 


where h 2 (x) and h 3 (x) are reduced functions and are given by the expressions 


h 2 (x) = ( qi -^-2q 2 s)ef^ dx 
ax 

h 3 {x) = (p- qi s + q 2 s 2 )e 2 ^^ dx 


( 14 ) 
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Using the lemma, the reduced Abel equation corresponds to a linear ODE without higher-order dissipative terms 

x^ + h 2 (x)x^ + h 3 (x) = 0. (15) 

On the other hand, in the case of f 3 = 0, the Abel equation of the first-kind ([3]) becomes a Riccati equation, while 
the second-kind homogeneous Abel equation 0 is reduced to the Riccati equation 


which is equivalent to 


^ = ~f 2 {x) - fi{x)v - f 0 (x)v 2 , 


^ = fo{x) + fi(x)y + f2(x)y 2 


(16) 


(17) 


This corresponds to the simple fact that an inverse power of a Riccati solution also satisfies a Riccati equation with 
redistributed coefficients. Now, we eliminate the linear part to obtain the reduced Riccati equation 

dz 


dx 


= h 0 (x ) + h 2 {x)z 


which corresponds to a nonlinear ODE with higher-order dissipative terms 

x ££ + h 2 (x)x<; + h 0 {x)x = 0 , 

where now the coefficients are 

h 0 (x) = re~^ q2 ~ 3rs)dx 

h 2 {x) = { qi -s'- 2 q 2 s + 3 rs 2 )e^ q2 - 3rs)dx . 


(18) 

(19) 

( 20 ) 


III. SOME ABEL INTEGRABILITY CASES 


A. Abel’s equation of constant coefficients 


Denote the coefficients fi{x) = Ai of m, where Ai £ R. are constants and A 3 ^ 0 so that F{x 1 y) = F{y). It is 
obvious that the roots of the equation F(y) =0 are themselves solutions of ([3]). More generally, the general solution 
of ([3]) is obtained via factorization of the denominator in the right-hand side of 


dy 


Aa y + A , 2 / + 


= A 3 dx , 


which leads to the following cases: 

(i) If yi ^ y 2 ^ y 3 £ R , so integration of (12T1) has the form 


( 21 ) 


dy 


J (2/-2/1X2/-2/2X2/-2/3) 

(ii) If yi ^ 2/2 = 2/3 £ R, then the integration of (HHI) is of the form 


= (y-yi ) V2 V 3 (y~y2 ) V3 V1 (y - 2/3X 1 y2 =cie 


A 3 x 


( 22 ) 


dy 


1 ,i y-yi\ A 

In |-1 = A 3 x + c 2 


J (2/- 2/1X2/- 2/2X (2/1 - 2/2)(2/ - 2/2) (2/1 -i/a ) 2 22-2/2 

(iii) If 2/1 = 2/2 = 2/3 € R , so (l2Tjl becomes 


(23) 


dy 


J (2/ — 2/i ) 3 —2(2/ — 2/1) 2 

(iv) If 2/i = 2/2 = a + i/3 £ C and y 3 £ R, the integration of (12T1) leads to 


= A 3 x + c 3 , 


(24) 


dy 


((2/ - a ) 2 + / 3 2 ) (2/ - 2/3) 


= In 


V(y- a ) 2 + X 2 


a-yi , y-a 
H--— arctan —-— = A 3 x + c\ 


(25) 
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B. Integrability based on the normal form of Abel’s equation 

If the following transformations as given in Kamke’s book [l3| 


y(x) = u(x)r](£(x)) - 
t(x) = j /3 w 2 dx 


/ 2 W 

3/3(x) 


ui(x) = e 




(26) 


are applied to equation (|5|). then one obtains Abel’s equation in normal form 


where the invariant I(x) is given by 


I{x) = 


% = ” 3 + /(I) 


fn + - — (A — fl * 2 + 
3 dx \ f 3 J 3/3 


2 fi 
277/ 


/3W 3 


(27) 


(28) 


Thus, we conclude that if /( x) = constant, then (1271) is integrable since it is separable. If one chooses relations between 
the functions /, such that the invariant is null and letting f 2 {x)/f^,(x) = n(x), then 


^ + ^/ 2 « 2 - fin + 3/o = 0 


(29) 


is a Riccati equation, which is always integrable because it is obtained from Abel’s equation in normal form, which is 
integrable. Thus, the solution to (1271) is 


1 


V = 


V^2f ’ 

and explicitly, the solution to ([3]) with null invariant I ( x ) is 

/ (. h-m) dx 


(30) 


y = 


Ic — 2 J fs 


A 

2 3/ 3 


( 31 ) 


C. Integrability of Abel equation with non-constant coefficients and fo = 0 

In this subsection, we will consider the original Abel equation ([5]) of the first-kind 

^ = h{x)y + f2{x)y 2 + f 3 (x)y 3 , (32) 

which corresponds to ([4]) without the cubic nonlinearity. 

Let us use the transformation 

P f h{x)dx 

y ( x ) = -— . (33) 

v{x) 

which is © with z = i that allows us to put (1371) into a differential form 

vdv + (P + Qv)dx = 0 , (34) 

where 

P(x) = f 3 (x)e 2 f fl(x)dx 
Q(x) = f 2 (x)ef h ^ dx . 


(35) 
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Thus, the original equation m is considerably reduced only by having /o = 0. 
For equation (1351) . if there exists an integrating factor 0 



H(x,v) = e~ k ^ +dQ{x)dx) , 

(36) 

where k is a constant, then 

(P + Qu)ndx + pfidp = 0 , 

(37) 

is integrable, provided that 

Qy, + (P + Q v )^ = , 

(38) 

which leads to Q = kP. Hence, the necessary condition for the integration of (1521) is 



kh(x)e !M^ = Hx) _ -4(A) = /,£-> d £ = fin 

(39) 

Using /x(a;, u) = e kl, e fc2 

f p{x)dx, we thgt there exists the potential ^(x,p) = c$, which satisfies 



dT 

dx 

dT 

dv 


= (1 + kv)Py, , 


= p/i. 


(40) 


Therefore, the solution to (l34l) is 


*(x,v) = (1 + kv)e- k ( v+k f pdx )=c*, 


(41) 


or in terms of the integrating factor 


^(a;, v) = (1 + kp)/j,(x, v) = c^. 

By substituting the condition (l39l) . the invariant (l28ll has the form 


T , n 2fc 3 j tl dx 
I{x) = — e J * . 


(42) 


(43) 


D. Canonical form of Abel’s equation and the integrating factor 

According to the book of Kamke 0 , for equations of the type (1521) for which there is no condition (1551) , one should 
change the variables according to 


y(x) = u(x)f}( C(x)) , u> = e-T /l(x)dx , 
C,{x) = j f 2 {x)oj{x)dx , 


which lead to the canonical form 


- ~2 


d( 


= v +g{ Ov 


(44) 


(45) 


where 


ff(C(®)) 


fs( X ) r [ h(x)dx 

h(x) 


(46) 


is the Appcll invariant. Then, the integrating factor /r can be used to formulate the following interesting result: 
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Lemma 2. Any Abel equation in the canonical form (l45l) is integrable as long as the invariant is constant, with 
g(Q = i, and solution given by (TUI) . 

Proof. In HMD, let u ( x ) = ’ w h ere dC t = Q(x), which leads to 


df) 

d( 


P(x)_ 3 , h{x)eif^ dx n 1„ 3 

yyr4v 3 = V 2 + ’ , , - V 3 = V 2 + T V 3 

Q(x) J2\x) k 


that is separable, with solution 


1 , 1 1 

k ln fj + k 


c 


1 

V 


which after simplification leads back to m- 


(47) 


(48) 


IV. VEIN’S ABEL EQUATION 


The following Abel equation: 

dy 


-2b 


dx bx + a 2 


y- 


3 (ax + b 2 ) 2 x 3 — 3abx — a 3 — b 3 


bx + a 2 


-y 


bx + a 2 


-y 


(49) 


where a, b £ R, is known from a paper of Vein to be integrable We first notice that if we define the following 
determinants, 


D, = 


x —a —b 
—b x —a 
—a —b x 


D* = 


x —b 
—b —a 


D 1 = \-b\, D = 


—a —b 
x —a 


then they allow us to write the coefficients of Vein’s Equation (1491) as /i = 2Di/D, f 2 (x) = —3D 2 /D, and fs(x) = 
Dz/D. Thus, (147)1) can be written in the compact form 


dy _ 2D\ 3 D 2 2 D 3 3 

dx~ D V D V + D V ' 

Next, using the transformation given by (1261) with the coefficients from (1491) . we have 

a 2 + bx 

U! = 


x 3 — 3 abx — a 3 — b 3 


(50) 


(51) 


and 


n , , 2 

di = h{x) “ 


Using these expressions, the invariant I(x ) becomes 


(52) 


d f2 _ fl /2 _i_ 2/2 2 

, . dx 3/3 3/3 ^ 27 7 i -U 

1 (x) = - 1 —^-— = -—- = — 1 . 


/3W 3 


l w 3 


(53) 


Therefore, the normal form of Abel’s Equation C5l) is 


(iH 3 - 1 - 

dS, 


Using ESI) and integrating (15S1) . we obtain 


. 1, (77 — 1) 2 , 

£(»?) - c = - In 2 - — arctan — (1 + 2ri) . 

6 I + 77 + 77 2 3 3 


(54) 


( 55 ) 
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But because £(x) = f uidx, we find 


f(x) — c = — In_ ( X ^ (a + 6) ) 2 _ 

6 x 2 + (a + b)x + a 2 -ab+b 2 


V3 


arctan 


■\/3 2x -(- a b 
3 a — b 


From the last two equations, one might infer that 


( 56 ) 


and consequently using (1551) would get 


r/(x) 


x + b 
a — b 


y{x) 


b 

(x — (a + 6)) (a — b) 


( 57 ) 


( 58 ) 


Unfortunately, because for (1571) . the argument of arctan is satisfied, but the argument of In is not satisfied, then we 
conclude that (1551) is not the solution of (1151) . 

To see how the implicit solutions (l55l) . (l57l) can be untangled and written explicitly, we will proceed as Vein did. 
Firstly, we show that there are actually three solutions of the type (1551) . which are the solutions of (1151) and can be 
put in the form 0: 


y i 


= (ax + b 2 ) — (bx + a 2 ) 


2/2 


= (ax + b 2 ) — (bx + a 2 ) 


t 2 


2/3 


= (ax + b 2 ) — (bx + a 2 ) 


, ( 59 ) 


with the t-functions defined cyclically as 
ti = ccl)i(s) + (j> 2 (s) 


h = C(f> 2 (s ) + (s) , t 3 = C(j) 3 (s) + 4>i(s) , 


where c is an arbitrary constant as it does not appear in the differential equation. The </> functions are the following 
triad: 

Mx) = \[ &X + 2e _a:/2 cos , 


M x ) [e x ~ 2e x/2 cos (^- + ^)] 
Mx ) [e x - 2 e ~ x/2 cos (^- - ^)] 


( 60 ) 


also known as the third-order hyperbolic functions, and are plotted in Fig. [TJ 



FIG. 1: The third-order hyperbolic functions <j>i, i = 1, 2, 3 of Equation @9. 
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They are independent because their Wronskian W(0i, 0 2 , 03) = 1, and they satisfy the following relationships 

f <£i + 02 + 03 = e x , 

1 01 + 02 + 03 - 3010203 = 1 • 

They also fulfill the following relationships 


d(j)i/dx = 03 , 

d 2 <pi/dx 2 = 02 , 

d 3 cpi/dx 3 = 0 i , 


dcf)2/dx = 0 i , 

d 2 (p 2 /dx 2 = 03 , 

d 3 (p 2 /dx 3 = 0 2 , 

(61) 

dcj )3 / dx = 02 , 

d 2 (p 3 /dx 2 = 0i , 

d 3 (p 3 /dx 3 = 03 , 



which are cyclic. That is, the fourth-order derivatives have the same property as the first ones, the fifth-order 
derivatives as the second ones, and so forth. In particular, they are independent solutions of the differential equation 
d 3 (p/dx 3 = 0, but also of d 6 (p/dx 6 = 0 and of any d 3+3m (p / dx 3+3m = 0, with m = 0,1, 2,.... 

Let 


t = C 0 i(s) + 02(s) , 

where c is a constant. Then, t fulfills 

t' = c0 3 (s) + 0i (s) , 
t" = C0 2 (s) + 03 (s) , 
f'" = t 

with ' = d/ds. Suppose s is defined implicitly as a multi-valued function of x as follows 


f 


x = a[ — +6 — = 


a0i(s) + 6 c0 2 (s) + (ac + b)cp 3 (s) 
C 0 i(s) + 02 (s) 


= $(s) 


and denote the inverse relation by the explicit formula 

s = $ -1 (x) . 

Then y = ds/dx = d^~ 1 (x)/dx satisfies Abel’s equation C^l) . 

Indeed, from (l64l) . differentiating with respect to s and by the usage of (l63l) . one obtains 


dx 

ds 


t 




The quotient L- can be eliminated by means of (RHl) . which leads to 

b— = (ax + b 2 ) — ( bx + a 2 )— . 
as t 

Differentiating again, this time with respect to x we have, 


dx 2 


(—V 

V dx ) 


= a — b 


— (bx + a 2 ) 


2 1 


ds 

dx 


(62) 

(63) 

(64) 

(65) 


( 66 ) 


Eliminating the t-quotients by means of (1641) and (1661) . one can find after some reduction the following equation 


dr s 


ds 


(bx + a 2 )—— = —2b— + 3(ax + b 2 ) ( — + (x 3 — 3 abx — a 3 — b 3 ) 


dx 2 


dx 


ds \ ‘ 


dx J 


’ — a 3 — b 3 ) (-j— 
\ dx 


Substituting V = 'jr and dividing by (bx+a 2 ) lead to Abel’s equation C9l) . Thus, the solution can be written explicitly 


as 


- = (ax + b 2 ) — (bx + a 2 ) 
V 


c 0 3 ( 4 » 1 (x))+ 0 i( 4 » 1 (x)) 
c0i($~ 1 (x)) + 0 2 ( 4 >~ 1 (x)) 


(67) 


with c an arbitrary constant. 

The other two solutions can be obtained by cyclically replacing the 0’s in the t function 

h = c0 2 (s) + 03 (s) , t 3 = c0 3 (s) + 0i (s) . 


( 68 ) 
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V. VEIN’S NONLINEAR OSCILLATOR 


To see what kind of oscillator corresponds to Vein’s Abel equation, we will proceed as follows. First, we will 
eliminate the linear term from (1491) using 


to obtain 


-2b f 

y — ze bx + a z 


z 

(i bx + a 2 ) 2 ’ 


(69) 


where 


dz 

dx 


h 2 (x)z 2 + h 3 (x)z 3 , 


h 2 (x) 

h 3 (x) 


3 (ax + b 2 ) 

(bx + a 2 ) 3 
x 3 — 3 abx — a 3 
(bx + a 2 ) 5 


b 3 


Then using lemma 1 with xj = v(x(£)) = we obtain 


(70) 


(71) 


x^ + h 2 (x)x£ + h 3 (x) = 0. (72) 

which represents a nonlinear oscillator with rational friction and nonlinearity that are plotted in Figures [2] and (3( 
respectively, for the following values (i) a = 1, b = —2; (ii) a = 1, b = —1; (iii) a = 2, b = —1; (iv) a = —1, b = 1. 


h 2 (x) 



(- 1 , 0 ) 

( 0 , 0 ) 

( 1 , 0 ) 

— ( 0 , 0 ) 


FIG. 2: The friction function h 2 (x) of Vein’s oscillator at the fixed points (—1,0), (0,0), (1,0), and (0,0). 

Writing equation (1721) as a dynamical system 

xq = v = M(x,v) , , 

= — h 2 (x)v — h 3 (x) = N(x, v) ^ ' 

and because ((731) cannot be put in the form of (l40l) . because ^ const., we will be using instead the standard 
methods of phase-plane analysis and use the linear approximation at the equilibrium points of (1731) to classify them 
instead of solving the equation by finding the potential 41. 

Also, notice that the system is non-Hamiltonian because there is no potential 't , (x,u) = const, such that 

If- ~ v 

H = h 2 (x)v + h 3 (x). 


(74) 
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The Jacobian matrix of (1731) is 


J = 

■ dM_ dM_ ' 
dx dv 


dN 

dN 



. dx 

dv . 



dh 3 

dx 


1 

- h 2 (x ) 


(75) 


We have three equilibrium points of the system (1731) from which one is real (xo,?/o) = (a + b, 0), while the other two 

are a pair of complex conjugates, (xi } 2 ,yo) = {cx±i( 3,0), with a = and /3 = — . The particular case a = b 

which Vein considered, will reduce the three fixed points to the case of one real fixed point (2a, 0). The characteristic 
polynomial of the Jacobian matrix is 


\J-\I 2 \ =A 2 -<5iA + J 2 = 0 , (76) 

while the discriminant is A = — 4^2- By evaluating the Jacobian at the real fixed point J = J\(c,o)i where C is 

either real or complex, we obtain 


from which 


J = 


0 1 

-h 2 (C) _ ■ 


(77) 


\J-\h\ 


X 2 + h 2 {C)X + ^{C) = 0 . 
dx 


1. Real fixed point: In this case we have 


(78) 


— (a*+ab+b^y J < 0 
^ 2r = (a^+ab+b^) 4 > 0 ’ 

A r = —$2r < 0 . 


(79) 


then, the real fixed point (a + 6, 0) living on the x axis is a stable spiral. The eigenvalues of the Jacobian matrix are 
complex conjugate pairs Ai ,2 = ^(5i ±iy/$ 2) = 2(a^+ab+b a ) 2 ± i). By calculating the eigenvectors, the linearized 

solution around the fixed point is 


z(C) 

HO 


= g 2(a 2 +ob+b^)^ 


Cl C0S 2 (a'l+ab+b*)* + C2 Sin 2 (a'l+ab+b'^ 

_ ~ (a' 2 +ab+b' 2 ) 2 { Cl S ^ n (2 (a' 2 +ab+b 2 )* + 3) _ C<1 C0S ( 2(a 2 +ab+b 2 ) 2 ~ 3 ) } . 


(80) 
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Fixed Points 

<5i 

82 

A 

Type 

{x 0 ,yo) = (a + b, 0) 

&lr 

&2r 

— S2r 

stable spiral 

{xi, 2 ,yo) = (--|[(a + b) ± W 3(a - 6)],0) 

Sic 

82c 

— 82c 

no conclusion 


TABLE I: Classification of equilibrium points of m- 


where ci and C 2 are arbitrary constants. 

2. Complex conjugated fixed points: In this case the coefficients of the characteristic polynomial <5!, 82 are complex, 
and their values are shown in Table U and given in (1811) . 

8i c = 2(a^-b‘ i )' i [( q2 - 2ab ~ 2 ^ 2 ) + iV3a(a + b )] , 

82c = 2 ( a 3 _fc 3 ) 4 5 6 7 8 9 10 [(~ ffl4 — 8 a 3 b — 6 a 2 b 2 + 4afo 3 + 26 4 ) + iy/3a(a 3 — 6 ab 2 — 4fe 3 )] , ( 81 ) 

A c = 82c • 

Because both coefficients 8 \ c and 82 c are not real nothing can be said about these complex fixed points. 

The phase-plane portrait of the system (l73l) is shown in Fig. [4] The real fixed point is identified by the red dot, 
and the isoclines by the dotted curves. 


VI. CONCLUSION 

This paper recalls several integrability properties of Abel equations together with some simple consequences, followed 
by a discussion of Vein’s Abel equation related to third-order hyperbolic functions. The corresponding nonlinear 
oscillator system is introduced here, and its dynamical systems analysis is provided. Finally, it is worth noticing that 
the nonlinear oscillators corresponding to Abel equations are in fact a large class of generalized Lienard equations of 
the form x + f{x)q{x) + h(x) = 0, where g , and h are arbitrary functions and f(x) = a\X 3 + a 2 X 2 + a 3 X + aq, where 
cti are constants |15j| and, as such, have many applications in physics, biology, and engineering GM3- 
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